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Abstract 

This is the second of two papers devoted to the perturbative computation of the 
ghost and gluon propagators in SU (3) Lattice Gauge Theory. Such a computation 
should enable a comparison with results from lattice simulations in order to reveal the 
genuinely non-perturbative content of the latter. The gluon propagator is computed 
by means of Numerical Stochastic Perturbation Theory: results range from two up 
to four loops, depending on the different lattice sizes. The non-logarithmic constants 
for one, two and three loops are extrapolated to the lattice spacing a — > continuum 
and infinite volume V — > oo limits. 

1 Introduction 

In [1] a three-loop computation of the ghost propagator was presented. As announced, we 
now report on a similar computation of the gluon propagator. The full propagators of glu- 
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ons and ghosts calculated on the lattice encode information about the non-perturbative vac- 
uum properties of QCD and of pure Yang-Mills theory [2]. This requires a non-perturbative 
extension of the Landau gauge. 

Although BRST invariance is essential for other non-perturbative approaches to non- 
Abelian gauge theory, there are principal difficulties to reconcile it with present-day lattice 
gauge fixing technology (as applied for calculating gauge-variant objects like gluon and 
ghost propagators) in a way avoiding the Neuberger problem (see [3] and references therein). 
So far, there is no generally accepted way to deal with the Gribov ambiguity in lattice 
simulations. 

Other non-perturbative methods for calculations in the Landau gauge like Schwinger- 
Dyson (SD) equations [I] or the Functional Renormalization Group (FRG) [5j E] are not 
explicitly taking into account (and seem not to be affected by) the complication compared 
to the standard Faddeev-Popov procedure that the functional integration should be re- 
stricted to the Gribov horizon. Zwanziger has argued [7J that the form of the SD equations 
will not be affected. Rather supplementary conditions restricting solutions should account 
for it. The so-called scaling solution [SJ E] has been an early result of the SD approach 
obtained for power-like solutions. The power behavior has been further specified in [10] 
and confirmed by (time- independent) stochastic quantization [TTj . 

Concerning lattice simulations, Zwanziger proposed [12J to add nonlocal terms to the 
action that should control the simulation to stay within the Gribov horizon Q. In the 
context of the Schwinger- Dyson approach this possibility has been analytically discussed 
first in Ref. [13]. Standard Monte Carlo simulations without these refinements have failed 
to reproduce the theoretically preferred far-infrared asymptotics of the scaling solution 
and have supported instead the so-called decoupling solution (see Ref. [H] and references 
therein). This solution has later been shown to be possible in the SD and FRG approaches 
with suitable boundary conditions [13], but at the expense of a conflict with global BRST 
invariance. 

Although we can assume that NSPT remains in the vicinity of the trivial vacuum, we 
have to understand the present situation with respect to Monte Carlo results for gluon 
and ghost propagators. Fortunately, the momentum range we are interested in and where 
we are going to compare with Monte Carlo lattice results is not influenced by the Gribov 
ambiguity and the way of gauge fixing [TB]. From various studies it is known, however, 
that the intermediate momentum range of (9(1 GeV) is not less important from the point 
of view of confinement physics, as seen for the gluon propagator [17J and the quark prop- 
agator [18] in Landau gauge on the lattice. From the FRG approach it is known that for 
the onset of confinement at finite temperature the mid-momentum region of the propa- 
gators is important [TJJJ |2U]. This is the region where violation of positivity [I7J EB l2"2"] 
invalidates a conventional, particle-like interpretation of the gluon propagator. Specific 
non-perturbative configurations (center vortices) have been found to be essential JT7J [18] 
to understand the behavior of the gluon and the quark propagator in the intermediate 
momentum range. 

In recent years, also the large-momentum behavior of the lattice gluon and ghost prop- 
^n this first attempt a singular ghost propagator oc 1/p 4 has been obtained. 
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agators has attracted growing interest, in particular with respect to the coupling constant 
of the ghost-gluon vertex J23J 1211 [251 12S] , which has the potential to provide an indepen- 
dent precision measurement of a s [Mz) from these propagators [27]. First estimates of 
the zero- and two-flavor values of [2H [251 EE] and a possible dimension-two conden- 
sate [211 [26] are available already and look promising. Therefore, a detailed knowledge of 
the propagators' lattice perturbative part would much foster these efforts. 

When our work begun in 2007 [28l [29] , the intention was to clarify the perturbative 
background, among other facts the type of convergence of the summed-up few-loop per- 
turbative contributions to the propagators in various momentum ranges. In standard 
Lattice Perturbation Theory (LPT) such calculations are very difficult beyond two-loop 
order. To overcome this obstacle, we have used the method of Numerical Stochastic Per- 
turbation Theory (NSPT) [30], which provides a stochastic, automatized framework for 
gauge-invariant and gauge-non-invariant calculations. 

Stochastic gauge fixing is built in, and a high-precision procedure has been devised 
to fix the gauge to Landau, at any given order. Thus, the propagators (for each lattice 
momentum) can be calculated in subsequent orders of perturbation theory. A limit is set 
only by storage limitations and machine precision. There are no truncation errors. For 
direct comparison with Monte Carlo results, we present the low-loop results summed up 
with inverse powers of the bare inverse coupling /3 n . In this paper we also try, for the first 
time applied to the gluon propagator, to improve the convergence by applying boosted 
perturbation theory. 

The effectiveness of NSPT relies on the fact that the parametrization of the (leading 
and non-leading) logarithmic terms can follow largely in accordance with standard pertur- 
bation theory. The essential difficulty left to NSPT is the computation of the constant 
contributions, which are in general very difficult to achieve in diagrammatic LPT. 

Only the one-loop constant term was known since long [31] for the ghost and gluon 
propagator. Reproducing these results, which are obtained in the continuum and infinite- 
volume limits, was the first feasibility test for NSPT [28J. 

In general, at any order, NSPT results are obtained at finite lattice spacing and finite 
volume. A fitting procedure is needed to get the continuum (a — > 0) and infinite-volume 
(V oo) limits. While the extraction of the first (continuum) limit relies on hypercubic- 
invariant Taylor series [52], a careful extraction of the second (infinite- volume) limit re- 
quires the accounting of pL contributions (p being the momentum scale relevant to the 
computation and L the finite extent of the lattice). In the first paper (abbreviated as I), 
we gave a quite comprehensive description of all this technology while applying it to the 
ghost propagator. In the present, second paper we are going to apply the method to an 
analysis of the gluon propagator. 

The paper is organized as follows. Sect. [2] recalls the lattice definition of the gluon 
propagator, together with specific features that the calculation in the framework of NSPT 
contains. Sect. [3J contains the nomenclature of standard Lattice Perturbation Theory where 
our results have to fit in. In Sect. H]we only briefly describe the implementation of NSPT. 
The interested reader will find more information of this kind in part I of this series of 
papers. However, we document the statistics for different lattice volumes and different 
orders of perturbation theory that has been collected by the Leipzig and Parma part of 
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the collaboration. Here we also present the raw data before and after the extrapolation 
to the Langevin time-step e — > limit. In Sect. [5] we compare the results of NSPT (up 
to four loops) with Monte Carlo results and try to improve the convergence by boosted 
perturbation theory. Sect. E] presents the fitting procedure and the final results for the 
leading and non-leading loop corrections. In Sect. [7] we draw our conclusions and summarize 
our results. 



2 The lattice gluon propagator 

The lattice gluon propagator D ab (p(k)) is the Fourier transform of the gluon two-point 



function, i.e. the expectation value 

DfMk)) = (A;(k)A b u (-k)) = 5 ab D,Mk)) , (i) 

which is required to be color- diagonal and symmetric in the Lorentz indices fi, v. For 
the definition of the lattice momenta p^(k^), k^ and p^(k^) to be used later we refer to 
(I-17)-(I-19). 

Assuming reality of the color components of the vector potential and rotational invari- 
ance of the two-point function, the continuum gluon propagator has the following general 
tensor structure § 

DM - - >-f) o(r) + , (2) 

with D(p 2 ) and F(p 2 ) being the transverse and longitudinal propagator, respectively. The 
longitudinal propagator F(p 2 ) vanishes in the Landau gauge. 

The lattice gluon propagator D^ v {j>{k)) depends on the lattice four-momentum p(k). 
Due to the lower symmetry of the hypercubic group its general tensor structure can be 
expected to be more complicated than (T5]) that holds in the continuum. Inspired by the 
continuum form fl2]) we consider as one strategy only the extraction of the following lattice 
scalars 

that should survive the continuum limit. Note, however, that additional lattice scalars 
could be measured as well. The first scalar vanishes exactly in lattice Landau gauge. In 
this gauge the second scalar function, corresponding to the transverse part of the gluon 
propagator in the continuum limit, is denoted by 

D(p(k)) = \i2 D »M k ))- (3) 

On the lattice, this function is influenced by the the lower symmetry of the hypercubic 
group. 



2 Diffcrcntly to the other chapters p denotes here directly the continuum Euclidean four-momentum. 
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In NSPT the different loop orders n (even orders in j3 1 ^ 2 ) at finite Langevin step size 
e are constructed directly from the Fourier transformed perturbative gauge fields A^ l '(k) 
with (see (1-6) and (1-7)) 

l>0 a 

what leads to 

/ 2n+1 „ \ 
8 ab D$(p(k)) = { J2 [ A f\k) Al^ +2 ~ l X-k)} j . (5) 

Note that already the tree-level contribution to the gluon propagator, D^J , arises from 
quantum fluctuations of the gauge fields with 1 = 1. In addition, terms with non- integer n = 
1/2, 3/2, . . . in the previous equation <^ - which do not correspond to loop contributions 
- should vanish numerically after averaging over configurations. 

Similar to the ghost propagator in paper I, we present the various orders of the gluon 
dressing function (or "form factor") in two forms: 

J G ^ n \p) = p 2 D {n \p(k)) , J G > {n \p) = p 2 D {n \p{k)) . (6) 

Contrary to the ghost propagator the gluon dressing function can be calculated at the 
same time for all possible four-momenta given by four integers - the four-momentum tuples 
(ki, &2, k%, ki). This makes a calculation of the gluon propagator significantly cheaper. The 
tree-level result for the dressing function, J G '^(p(k)) = 1 in the limit £ — > for all sets of 
tuples is non-trivial and is obtained as the result of averaging. 



3 The propagator and standard Lattice Perturbation 
Theory 

As discussed in paper I - we relate infrared singularities encountered in our finite volume 
NSPT calculation to powers of logarithms of the external momentum obtained in the 
infinite volume limit. Therefore, we need the anomalous dimension of the gluon field 
Ap, the fl- function and the relation between lattice bare and renormalized coupling. The 
procedure is outlined in detail in Section 4 of I. Here, we only repeat the essential equations 
and quote the final numbers. To avoid a possible mismatch of equations we add an index 
G for the gluon propagator case. 

In the RI'-MOM scheme, the renormalized gluon dressing function J G > RI is defined as 

tG,ri' /„ „ \ J G (a,P, am') m 



with the standard condition 



J G ' RI '(p,/i,a R i/)| p 2 = ^2 = 1 . 



The gluon dressing function J G (a,p, ccrf) is the gluon wave function renormalization con- 
stant Z G (a, fi, Ori') at jj 2 =p 2 . 
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The expansions of Z G ', J G and J G ' RI ' in terms of the renormalized coupling a(a/j) = 
(g(a/i)) 2 /(167r 2 ) are completely analogous to (I-39)-(I-41). The gluon wave function renor- 
malization as an expansion in the lattice coupling «o = N c / (87r 2 /3) is represented as 



(a. 



11, a ) = 1 + E «o E Z Z ( lo gM) J • (9) 

j>0 j=0 



This is the expansion we can measure in NSPT. 

Again we restrict ourselves to three-loop expressions for the Landau gauge in the 
quenched approximation. The coefficients in front of the logarithms are partly known 
from calculations of the gluon wave functions and the beta function in the continuum and 
given as follows (compare e.g. [33J) 

z?f' = -fN c , (10) 

z&r—y^ 2 , ^ RI ' = -^AT c 2 + 3^o RI 'iV c , (11) 



G.RI' _ N 3 G,RI' _ _5495 3 31 G,Rl' ^2 n ^ 

' 3 ~~ 18 ' ' 2 ~~ 36 2 '° ' 

<7rt' 2127823 „ T o 361 ... „ T a 1279 r?RT' »t9 31 r?RT' » T 
= AT 3 -| (3 AT r 3 Ar H 4o N c . 

3,1 3g88 c T g SL J c 10g 1,0 c -r 3 2,0 

The finite one-loop constant z 10 is known from the gluon self energy in standard infinite 
volume LPT [31] with the value 

z G RI ' = 60.3673076 . (13) 



Again the constants z G ^ 1 and z G Rl were not known so far. The form of the coefficients 
of Z G (a, fi, ao) in ([9]) up to three loops can be read off directly from equations (I-51)-(I-53) 
replacing there all 2y by z G y 

As a result, we present the gluon dressing function as function of the inverse lattice 
coupling (3 to that order 

J G < 3 - lo °P(a,p,/3) = 1 + E ^ G ' (i) (a,p) , J G ^(a, P ) = E Jg (log(ap) 2 ) J , (14) 

i=l " j=0 

with the one-loop coefficients 

J G X = -0.24697038, 

J G Q = 0.03799544 z G Q R1 ' = 2.29368, (15) 

the two-loop coefficients 

J 2 G = 0.08210781, 

J 2 G = -0.917978574 - 0.00938375 z G m ' = -1.484450 , (16) 
J 2 G = 0.10673710 z G ' Rl ' + 0.00144365 z G ' R1 ' = 6.443431 + 0.00144365 z G ' RV , 



z 



and the three-loop coefficients 



j£ 3 = -0.02963736, 

J 3 G 2 = 0.62856608 + 0.00311972 zff V = 0.81689534 , 

J 3 G X = -4.06861256 - 0.06123991 zff 1 ' - 0.00035654 zgf 1 ' 

= -7.7655009 - 0.00035654 zg, M ' = -6.174164 - 0.246970 J 2 G , (17) 

J 3 G = 0.375990 zff 1 ' + 0.00811105 z^f 1 ' + 0.0000548523 zff 1 ' 

= 22.697532 + 0.00811105 z^f 1 ' + 0.0000548523 zff 1 ' . 

Let us repeat it once more: the leading logarithmic coefficients for a given order can be 
exclusively taken from continuum perturbative calculations. The non-leading log coeffi- 
cients are influenced, however, by the finite lattice constants from corresponding lower 
loop orders. 

4 Results of the NSPT calculations 

4.1 Statistics 

To obtain infinite volume perturbative loop results at vanishing lattice spacing, we have 
to study again the limit e — > and different lattice sizes N. We have used N = 6, 8, 10, 12 
and N = 16, 20, 32 and studied the maximal loop order for the propagator n max = 4 and 
n max = 2, respectively. The accumulated statistics for the different e's and lattice sizes are 
collected in Tables [1] and [2j 



e 


N = 6 


N = 8 


iV = 10 


N = 12 


N= 16 


0.01 


1500 


750 


3000 


1000 


3000 


0.02 


1000 


750 


2000 


1000 


2000 


0.03 


1000 


750 


2000 


1000 


2000 


0.05 


1000 


750 


2000 


1000 


2000 


0.07 


1000 


750 


2000 


1000 


2000 



Table 1: Number of gluon propagator measurements up to four loops (N < 12) and up to 
one loop (N = 16) using the Leipzig NSPT code. 

The Landau gauge was defined by the condition (1-14). In the gluon propagator case we 
have used in the gauge fixing condition (1-20) Z max = 10 orders of the perturbative gauge 
fields to obtain the propagator up to four loops in the case of the smaller volumes, and 
used Z max = 6 for the bigger volumes with iV = 16, 20, 32. 

4.2 Raw data and check of vanishing contributions 

In Fig.CQwe present as an example the measured dressing function J G '(") as function of a 2 p 2 
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e 


iV = 16 


iV = 20 


iV = 32 


0.010 


7436 


5965 


810 


0.015 


3053 


3896 




0.020 


4725 


3015 


715 


0.040 


2827 


2633 


835 



Table 2: Number of tree-level, one- and two-loop gluon propagator measurements at lattices 
sizes iV = 16, 20, 32 using the Parma NSPT code. 
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Figure 1: Raw data J G ^ vs. a 2 p 2 at = 10 and e = 0.01. Left: Increasing loop 
contributions from below. Right: Vanishing non-loop contributions. 



for n = 0,1,2,3,4 and n = 1/2,3/2,5/2,7/2, respectively, for a lattice size 10 4 and the 
Langevin step e = 0.01. Note that due to the hypercubic symmetry different realizations 
of the four-momentum tuples k can lead to the same a 2 p 2 forming different branches. 

In the left part of Fig. [TJ the different contributions are shown starting from tree- 
level up to four loops. The higher loop contributions are of the same sign leading to a 
monotonic increase of the perturbative dressing function including higher and higher loop 
contributions. The right part of Fig. [1] nicely demonstrates the vanishing of the non- 
loop contributions. The size of the "approximate zeros" has to be compared with the 
corresponding sizes of the loop contributions of the "adjacent" integer loop orders. 

Fig. [2] demonstrates for a chosen n — ^ that the corresponding averages equal to zero are 
indeed reached at finite Langevin step size e. Typically, the errors increase with decreasing 
step size e and a 2 p 2 . 

In Fig. [3] we present the result of the linear extrapolation with e — > for the tree- level 
J G -(°) and one-loop J G '^ gluon dressing functions on a 16 4 lattice, together with the data 
of all inequivalent four-momentum tuples for e = 0.07 and e = 0.03, as a function of a 2 p 2 . 
Note to what precision the expected result "One" for J G '^ is reproduced. 

Figs. H]and|5] show how the one-loop and two-loop dressing functions J G ^ and J G>(n ) 
depend on a 2 p 2 and (ap) 2 , respectively. The different branches for the inequivalent off- 
diagonal four-momentum tuples in use are clearly seen, which are strongly volume depen- 
dent. 
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Figure 2: j G >W 2 ) vs. a 2 p 2 at iV = 16 and e = 0.01 and 0.07. 
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Figure 3: The e — >■ extrapolation of the tree level J G '(°) (left) and one-loop J G '^ 1 - ) (right) 
dressing function vs. a 2 p 2 at N — 16. Note that the one-loop result is shown in a smaller 
a 2 p 2 window. The e — > extrapolation is reached from above. 



5 The perturbative gluon propagator summed to four 
loops 

5.1 Naive summation 

Having obtained the different loop contributions to the dressing function in Landau gauge, 
we can sum up those contributions for a given inverse gauge coupling (3 to get an estimate 
of the perturbative gluon propagator. 

We calculate the perturbative gluon dressing function at a given lattice volume summed 
up to loop order n max for a given lattice coupling (3 as follows: 

J G (rw)= £ l^J G ' {n) (18) 

n=0 " 
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Figure 4: The one-loop dressing functions j G -^ (left) and J G '^ (right) at selected different 
volumes vs. small momentum squared a 2 p 2 and (ap) 2 - 
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Figure 5: Same as in Fig. H]for the two-loop dressing functions j G '^ and J G,<y2 \ 

and similar for J G (n max ). In contrast to the ghost propagator, here also the tree level 
"One" is calculated numerically (not shown in the Figures below) and we take the extracted 
numbers being near to the exact number one. 

Restricting ourselves to four-momentum tuples near the diagonal we obtain the fol- 
lowing results summarized in Figs. E] and It is interesting to notice that all loop 
contributions are of the same sign. Therefore, the summed dressing function most prob- 
ably represents a lower bound on the complete perturbative dressing function. We have 
to stress that our perturbative results at large a 2 ft 2 or (ap) 2 , being at the edges of the 
Brillouin zones, have nothing to do with the continuum limit in the ultraviolet but could 
eventually describe the Monte Carlo data in that large lattice momentum region. 
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Figure 6: The cumulatively summed perturbative dressing function J G (n max ) (left) and 
J (jimax) (right) up to four loops (two loops) using (3 = 6 at N = 8, 10, 12 (16, 20, 32). 
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Figure 7: Same as in Fig. [6] at (3 = 9. 

5.2 Summation using boosted perturbation theory and compar- 
ison with Monte Carlo results 

It is well-known that the bare lattice coupling g 2 = 6/(3 is a bad expansion parameter. 
Due to lattice artifacts like non-vanishing tadpole graphs, the lowest order perturbative 
coefficients are typically very large. Therefore, it has been proposed |34j to use a boosted 
coupling g\ instead of the bare g 2 to improve the relative convergence behavior. We use 
here a variant of boosting where g\ is computed via the perturbative plaquette P pert (g 2 ) 
calculated in NSPT (instead of using plaquette measurements from Monte Carlo simula- 
tions): 

d= p£m- (19) 



n 



Let us expand the perturbative dressing function J G (n max 

^max 

n=0 

The analogous expansion in the boosted coupling g 2 has the form 



in g 2 and write 



Ji 



G, 



2n :-G,(n) 



Inserting 



p 



pert 



n=0 



1+ T,9 2n c n 



(20) 



(21) 



(221 



n=l 



into ( fT9]) we compute the boosted expansion coefficients j G '^ as functions of the naive 
coefficients and c„. For the orders under consideration we thus obtain 



J 

:-g,(o) 
Jb 



Jb 



G,(l) 



Jb 



G,(2) 



:-o,(o) 

]G,(2) +CiJ .G,(l) 



(23) 



G,(3) 



Jb 



G,(4) 



•G,(3) 
•G,(4) 



+ 2 Cl j G ' (2) + (c? + c 2 )j G '« 



+ 3ci jG,(3) + (3c 2 + 2c2)j . 



G,(2) 



(cl + 3 ci c 2 + c 3 ) j 



G,(l) 



Of course, for the infinite series both expansions should coincide J G (oo) = J G (oo) - 
but for every finite, truncated order (n max < oo) the perturbative series differ. Since the 
plaquette is less than one, it is clear from f lT^j) that g 2 > g 2 . However, the boosted expansion 
coefficients become significantly smaller than their naive counterparts, I ^ \j G '^\- 

The combination of both effects results in an improved convergence behavior of the series 

(EU). 

An example of such a reordering using ( 123]) is given in Table |3]for J (n max = 4) assigned 



loop order n 


jG,(n) 


"•G,(n) 
Jb 





0.9804(85) 


0.9804(85) 


1 


0.3638(44) 


0.3638(44) 


2 


0.1995(31) 


0.07762(165) 


3 


0.1296(24) 


0.02445(61) 


4 


0.09070(193) 


0.007847(323) 



Table 3: Expansion coefficients j G '( n ) and j b defined in ( 1201 and (12 ip . respectively, for 
J G (n maj; = 4) assigned to the momentum tuple (1, 1, 1, 1) and the lattice size iV = 12. 

to the momentum tuple (1, 1, 1, 1) at lattice size iV = 12. For the expansion coefficients of 
the plaquette in NSPT we use the measured c n [35] 

Cl = -0.334998 , c 2 = -0.0337441 c 3 = -0.0137452 c 4 = -0.00729851 (24) 
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for the same lattice size N = 12 (neglecting their tiny errors). The two-, three- and 
four-loop expansion coefficients given in Table |3] are significantly smaller using boosted 
perturbation theory. 

In Figs. E] we show the effect of boosting for the dressing function J G (n max = 4) at 
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Figure 8: Comparison of naive and boosted perturbation theory for the gluon dressing 
function at iV = 12 with Monte Carlo data including real and complex Polyakov loop 
sectors (non-trivial) and restricted to the real (trivial) Polyakov sector from [36] at two 
different (j values. To show the non-perturbative low momentum behavior, also iV = 32 
Monte Carlo data are included at — 6. 

two different (3 values and compare it with Monte Carlo results of the Berlin Humboldt- 
University group [36] which has used the same gauge field definition (1-6) as implied in 
NSPT. The boosted result (n max = 4) is shifted towards the Monte Carlo results. Note 
that at the larger j3 value the Monte Carlo results significantly depend on whether the 
measurements are restricted to the real (trivial) Polyakov sector. We observe a reasonable 
agreement of the NSPT summed dressing function at the largest (a 2 p 2 ) with the Monte 
Carlo data obtained in the trivial Polyakov sector for all four directions (if necessary, 
reached by Z(3) flips). For details, see [36] . 

In contrast to Monte Carlo results [T6], we do not see even in higher- loop lattice per- 
turbation theory any sign for a suppression for J in the infrared direction (at small a 2 p 2 ), 
such that this principal difference must be related to non-perturbative effects responsible 
for the expected confinement behavior of the gluon propagator. One of interesting ques- 
tions might be to what extent the difference could be effected by some phenomenological 
constants. The gluon condensate might be a prominent example. 

6 Fit results for LPT in the infinite volume limit 

Analogously to paper I we extract the finite constants Jf for loop order i in the pertur- 
bative expansion of the lattice gluon dressing function f lT4|) in Landau gauge. We use the 
same strategy as outlined there. 
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First we subtract all logarithmic pieces (supposed to be universal and known) from the 
gluon dressing function for each momentum tuple and for all lattice sizes. Next we select 
a range in (ap) 2 = Y,(aP/i) 2 , a Pn = k^(2-n/N) with (ap) 2 min < (ap) 2 < (ap) 2 m&x . Within 
that range we identify a set of four-momentum tuples (k±, A; 2 , k 2 , k^) which is common to 
all chosen lattice sizes. The data in that set are assumed to have the same pL = p (aN) 
effects for each given momentum tuple. Since finite- volume effects decrease with increasing 
momentum squared, we choose as reference fitting point - for an assumed behavior at 
N = oo - an additional data point at (ap) 2 ~ (ap)^ ax from the largest lattice size at our 
disposal. Then we perform a non-linear fit using all data points from different lattice sizes 
N in that set plus the reference point correcting for finite size (no functional form guessed) 
and assuming a specific functional behavior for the H(4) dependence. This functional form 
analogously to (1-62) is a hyper cubic-invariant Taylor series (with (ap) n = X^Q'iV)"') 

j9 Q (ap) = Jg, + Ci,i (ap) 2 + c ii2 + c li3 (ap) 4 + c iA {(ap) 2 ) 2 + c i>5 + • ■ • . (25) 

Here J^ (ap) denotes the part of the dressing function which does not depend on pL = 
p (aN) and logarithmic effects. Finally we vary the momentum squared window and find 
an optimal x 2 region which allows us to find the "best" Jf Q . 

This time we use our NSPT data - extrapolated to zero Langevin step - from lattice 
sizes {N = 8, 10, 12, 16, 20, 32} for the one- and two-loop case. We take into account tuples 
with \ki — kj\ < A£; max similar to (1-58). In addition to the ghost propagator case some 
more tuples are allowed, e.g. A£; max = 2: (2,0,0,0), (3,3,1,1), . . . 

Since we have now a larger lattice volume available when setting the scale as explained 
in I, the possible optimal fitting window in the momentum squared (ap) 2 becomes larger 
and the number of four-momentum tuples fitted together from all lattices is increased. 
To get a rough estimate of the three- loop finite constant, we also use the available NSPT 
measurements from smaller lattice sizes {N = 6, 8, 10, 12}. 

As examples of how the fitting works we present in Figs. IHJ EH] and [TT] the handling 
of one-loop, two-loop and three-loop gluon dressing functions at a particular x 2 value 
switching off the estimated lattice artifacts in two steps: first the (pL)-effect then the H(4) 
dependence. 

The results for the one-loop constant for ten best x 2 -values using A/c max = 2 are shown 
in Fig. [T2] which agree within errors. The corresponding best two- and three-loop constants 
for AA; max = 2 are given in Fig. [131 Note that for individual x 2 values the momentum 
squared windows differ, depending on the number of contributing four-momentum tuples 
which enter into the non-linear fits. 

Table H] contains a summary of the non- logarithmic constants Jf Q we have found using 
three different selection criteria AA; max = 1,2,3 for the four-momentum tuples. The errors 
are estimated by equally weighting the mean deviations squared of both the individual fits 
and the sum of the "best" ten x 2 values. The results from the different criteria coincide 
within errors. We have obtained a very good agreement at the level of 0.5 percent with 
the expected exact one- loop result Jp given in f [T5|) . 

Comparing the accuracy reached here with that of the ghost propagator, let us mention 
again that in the gluon propagator case already the tree level is calculated from quantum 
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Figure 9: Fitting of the one- loop non- logarithmic coefficient. Left: Data points after 
subtracting logarithms for various lattice sizes compared to points using the non-linear fit 
of type (1-66) with ( 1251) analogous to (1-62). Right: Stars denote data after correction for 
finite- volume effect, the finite volume effect, i.e. they represent J^ (ap); full circles are fit 
points after correcting both finite volume and some hypercubic effects with the exception 
of those proportional to the coefficients c^i and in (1251) . 
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Figure 10: Same as in Fig. |9]for the two-loop non-logarithmic coefficient. 

fluctuations of the gauge fields. So the one-loop accuracy in the gluon dressing function 
case should be fairly compared to the two-loop accuracy of the ghost dressing function. 

To estimate the influence of the missing larger volumes for the three-loop constant, let 
us compare the two-loop constants obtained in the lattice volume sets (N = 8, . . . , 32) 
(one- and two-loop) and (N = 6, . . . , 12) (three-loop). The results are collected in Table 
From the numbers given there we conclude that the missing data sets in the three-loop 
fit of (J G '( 3 '°)} do not entail a significant change. There is a small tendency to somewhat 
larger numbers using larger volumes. This has to be taken into account as a systematic 
effect in our estimate of J^ . 
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Figure 11: Same as in Fig. |9]for the three-loop non-logarithmic coefficient. 
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Figure 12: One-loop results J± Q allowing Afc max = 2 for the smallest available \ 2 together 
with the exact value 2.29368 (see ( !T5|) ). 





Ak- 

""'max 


1 2 3 


(^i G o) 

(4: > 


2.318(56) 2.303(34) 2.292(69) 
7.939(123) 7.927(117) 7.897(112) 
31.61(33) 31.66(47) 31.56(46) 



Table 4: Non-logarithmic constants ( Jf ) of the gluon dressing function averaged over the 
best ten x 2 values of the non-linear fits using different selection criteria A/c max for the 
four-momentum tuples. 

Finally we decided to take the selection criterion Ak max = 2 as the most suitable one 
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Figure 13: Similar to Fig. [12] for the analytically unknown two-loop and three-loop non- 
logarithmic coefficients J^ and Jff . 



Ah 


<J?o> 


<J?o> 




N = 8, ... ,32 


N = 6,...,U 


1 


7.939(123) 


7.919(79) 


2 


7.927(117) 


7.876(106) 


3 


7.897(112) 


7.888(115) 



Table 5: Comparison of the two-loop constant ( J%, ) using two different lattice volume sets. 



and present our numerical results for the unknown non-logarithmic constants in the gluon 
dressing function of infinite volume lattice perturbation theory in Landau gauge: 



J 2 G = 7.93(12) 
J 3 G = 31.7(5). 



(26) 
(27) 



Collecting all results we can write (|14p in a numerical form (restricting to at most five 
digits after the decimal point) 



J G ' 3 - loop (a,p, 0) = 1 + - (-0.24697 log(ap) 2 + 2.29368) + 



(0.08211 (log(ap) 2 ) 2 - 1.48445 log(ap) 2 + 7.93(12)) + 



(2f 



1 / - ■ /- ' \2\ 3 , n 01 con (\ I \2\ 2 o i oi o > l - >2 



, v -0.02964 (log(ap) 2 ) + 0.81689 (log(ap) 2 ) - 8.13(3) log(ap) 2 + 31.7(5)) . 



A transformation to the RI' scheme can be performed using the relations given in Section [3j 
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7 Summary 



In the present work we have applied NSPT to calculate the Landau gauge gluon propagator 
in lattice perturbation theory up to four loops. The summed gluon dressing function is 
compared to recent Monte Carlo measurements of the Berlin Humboldt University group. 
Both (NSPT) perturbative and non-perturbative results are in terms of one and the same 
definition of the gauge fields, both in Landau gauge fixing and in measurements of the 
propagator. To improve the comparison, we have also summed our results in a boosted 
scheme showing better convergence properties. 

The key goal of the lattice study of propagators is to reveal their genuinely non- 
perturbative content, which asks for disentangling perturbative and non-perturbative con- 
tributions. The commonly used procedure goes through the fit of the high momentum 
tail by continuum-like formulae (anomalous dimensions and logarithms can be taken from 
continuum computation). While this lets us gain intuition, it opens the way to further 
ambiguities, since irrelevant effects give substantial contributions to the perturbative tail. 
At large lattice momenta our calculations indicate that the perturbative dressing function 
constructed by means of NSPT with more than four loops will match the Monte Carlo 
measurements, thus enabling a fair accounting of the perturbative tail. The strong differ- 
ence which is left over in the intermediate and - moreover - the infrared momentum region 
should then be attributed to non-perturbative effects. Power corrections [21] and contribu- 
tions from non-perturbative excitations [T7] 137] are serious candidates for the description 
of these (better disentangled) effects. 

The one-loop result for the perturbative gluon propagator of Lattice SU(3) in covariant 
gauges (and in particular Landau) has been known for a long time. Using our strategy 
for a careful analysis of finite volume and finite lattice size effects we find good agreement 
with this result. In ( )28|) we have summarized our (original) two- and three- loop results. 
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